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Abstract 

We consider the cosmological evolution induced by the free energy F of a gas of maximally su- 
persymmetric hctcrotic strings at finite temperature and weak coupling in dimension D > 4. 
We show that F, which plays the role of an effective potential, has minima associated to 
enhanced gauge symmetries, where all internal moduli can be attracted and dynamically 
stabilized. Using the fact that the heterotic/type 1 S-duality remains valid at finite tempera- 
ture and can be applied at each instant of a quasi-static evolution, we find in the dual type I 
cosmology that all internal NS-NS and RR moduli in the closed string sector and the Wilson 
lines in the open string sector can be stabilized. For the special case of D = 6, the internal 
volume modulus remains a fiat direction, while the dilaton is stabilized. An essential role is 
played by light D-string modes wrapping the internal manifold and whose contribution to 
the free energy cannot be omitted, even when the type I string is at weak coupling. As a 
result, the order of magnitude of the internal radii expectation values on the type I side is 
v^Ai a', where Aj is the ten-dimensional string coupling. The non-perturbative corrections to 
the type 1 free energy can alternatively be described as effects of "thermal El-instantons", 
whose worldsheets wrap the compact Euclidean time cycle. 
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1 Introduction 



The 5*0(32) heterotic and type I strings are dual perturbative descriptions of the same 
underlying theory This is easily observed at the level of the low energy effective actions, 
which are equivalent after certain field redefinitions. This follows from the fact that short 
massless supermultiplets have protected masses and that A/io = 1 supergravity coupled to ten 
dimensional super Yang-Mills theory with given gauge group is unique at the two- derivative 
level. One interesting facet of the equivalence is that in ten dimensions, the heterotic and type 
I string couplings are inverse to one another, and thus one has the opportunity to uncover 
strong coupling effects. In dimension -D > 7 (-D < 5), this leads to a strong-weak (weak- 
weak) duality, while ioi D = 6, string couplings and internal volumes are interchanged [2]. 

In the literature, most of the applications of string dualities have been based on BPS 
states and therefore restricted to models where supersymmetry is preserved in static uni- 
verses. In general, extending these ideas to non-supersymmetric cases (see jS] for some 
examples) and cosmological evolutions is difficult. However, such a project can still be ad- 
dressed within the context of no-scale models |4]. The latter are defined at the classical level 
by backgrounds associated to vanishing minima of a scalar potential, which admit a fiat di- 
rection parameterized by the scale of spontaneous supersymmetry breaking. The non-trivial 
vacuum energy, which arises at the quantum level, backreacts on the flat and originally static 
universe, and induces a quasi-static time evolution in the background fields [s]. 

To be specific, start with a dual pair of supersymmetric heterotic and type I models. As 
follows from the adiabatic argument of [6], one may implement on both sides a spontaneous 
breaking of supersymmetry, thus giving a new dual pair. For example, if the heterotic 
theory is in a perturbative regime and the spontaneous breaking at the classical level is 
compatible with flat Minkowski space, the cosmological evolution induced at the one-loop 
level can be reinterpreted in the dual type I regime. In this paper, we spontaneously break 
supersymmetry by considering the models at finite temperature. This can be implemented 
at the level of the two dimensional CFT by compactifying the Euclidean time on a circle, 
whose boundary conditions depend on the fermion number [?]. In this case, the one-loop 
heterotic effective potential discussed above is nothing but the free energy of a perfect gas 
of supersymmetric strings. Applying the heterotic/type I duality, we find the existence of 
novel contributions to the type I effective potential coming from light D-strings. Despite 
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being non-perturbative, these corrections have a large impact on the cosmological evolution, 
as well as on the low energy spectrum of the theory, even at weak type I string coupling. 

A second method to spontaneously break supersymmetry is by introducing "geometric" 
fluxes along internal cycles [s]. When the R-symmetry charge associated to the flux is the 
fermion number, this method is related to the finite temperature case by a double Wick 
rotation. In this paper, we only explore the thermal breaking for simplicity and clarity, as 
most of our results have a direct generalization to the second case. In realistic situations, 
one must include zero temperature spontaneous supersymmetry breaking before switching 
on finite temperature. In this case, a general picture arises, where the induced cosmology 
can be divided into different stages. In the Hagedorn era, where the temperature T is close 
to the string scale Ms, a phase transition between pre- and post-big bang evolutions takes 
place. It can be described along the lines of Refs [9 10 at the level of the two dimensional 
CFT and is both free of initial singularity and consistent with perturbation theory. As the 
temperature drops, the cosmology induced by the one-loop effective potential can be trusted 
until infrared effects become relevant, such as in the cases of radiative breaking or confining 
gauge groups. For example, in standard GUT scenarios, this defines intermediate eras where 
the temperature evolves in either of the ranges Ms > T > Agut ot Aqut > T > Mew? 



where Aqut and Mew are the GUT and electroweak scales 11-14 . These intermediate eras 
are connected by a phase transition where the dynamics responsible for the breaking of the 
GUT group must be taken account. After the electroweak phase transition, the conventional 
history of the universe follows with the hadronic, leptonic and nucleosynthesis eras... 

One feature of the above Hagedorn and intermediate eras is the possibility to stabilize 
internal moduli (T2|[l5|[T6]. This is an important issue since current observations of the 
gravitational force place lower limits on scalar masses (see for example (l^). Many ap- 
proaches address this question by considering compactification spaces where (geometrical or 
non-geometrical) internal fluxes are switched on at the outset, while preserving some amount 
of supersymmetry [is]. This leads to a partial stabilization since fiat directions always per- 
sist in such models, at least at the perturbative level. However, we would like to stress that 
once supersymmetry is broken, fiat directions are generically lifted in string theory. This 
was considered long ago in non-supersymmetric heterotic string backgrounds, such as the 



SO (16) X 5*0(16) tachyon free theory toroidally compactified 19 . However, minimization of 
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the moduli-dependent "cosmological constant" generated by loop corrections in such models 
leads to an unacceptably large vacuum energy at the minima, since supersymmetry is ex- 



plicitly broken at the string scale. In 20 , it was realized that a gas of string modes, which 
carry both winding and momenta, generate a free energy that enables stabilization of radii 
moduli. Upon introducing a zero temperature spontaneous breaking of supersymmetry at 



the string tree level, it was shown in 13-16,21 that this effect also has a quantum version, 
with the thermal gas and free energy replaced by virtual strings which induce an effective 
potentiaQ An advantage of this type of stabilization is that during the intermediate eras, 
the induced masses are not constant. Instead, they follow the time-evolution of the tem- 
perature T{t) and supersymmetry breaking modulus M{t), which drop proportionally. It 
is only after the electroweak phase transition that M{t) is stabilized and that the induced 
moduli masses become constant. As a result, the energy of the moduli with time-dependent 



masses is diluted during the intermediate eras, and the cosmological moduli problem [23|j 
is avoided. Moreover, the decrease of M(t) oc T(t) during the intermediate eras gives a 
dynamical explanation of the hierarchy between the supersymmetry breaking scale and the 
string scale, M -C M^. 

This above dynamical moduli stabilization relies on the existence of perturbative states 
in the string spectrum, whose masses are determined by the expectation value of the moduli 
and vanish at the stabilization points. For instance, in toroidal or orbifold compactifications 
of the heterotic string, if the radius Ri of some factorized internal circle is not participating in 
the spontaneous breaking of supersymmetry, it can be attracted to the self dual point Ri = I 
associated to an enhanced SU (2) level one Kac Moody algebra. Another simple example can 
be realized in type II superstring, when the internal circle is used to spontaneously break the 
supersymmetries generated by the right-moving sector via the Scherck Schwarz mechanism. 
In this case, Ri can be stabilized at the fermionic point Ri = l/\/2 corresponding to a 



In Refs 22 , the effect of the Coleman- Weinberg effective potential is explicitly subtracted in order to 
isolate the backreaction on the moduli arising from particle production near extra massless species points. 
To be substantial, this mechanism supposes the moduli already have non-trivial motions at tree level. Since 
the no-scale models are based on classical static backgrounds, the moduli velocities occur as backreactions 
of the one-loop effective potential and particle production is higher order in perturbation theory. 

simplified statement of this problem is that the energy of scalars with constant masses dilutes slower 
than the thermal energy of radiation, and so heavy scalars tend to dominate at late times, which can cause 
problems for nucleosynthesis. This may be fixed by requiring the heavy scalars to be unstable so that their 
fiuctuations eventually decay, thereby reheating space-time. However, the reheating process creates extra 
entropy and one can run into problems with baryogensis. 
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Kac Moody level two SU{2) extension [9]. However, since this type II setup is intrinsically 
left/right asymmetric, it cannot be extended to orientifold models in a straightforward way. 
Thus, the purpose of the present work is to infer how the internal moduli in type I no- 
scale models are stabilized by using our knowledge of the dual heterotic picture. As said 
before, we consider only thermal effects, as this is sufficient to uncover the mechanism. 
More specifically, using heterotic/type I duality at finite temperature, we infer the existence 
of non-pertrubative contributions to the thermal free energy of type I superstrings. These 
contributions are due to light, or even massless, D-strings which wrap the internal cycles 
and participate to the dynamical stabilization of all the internal moduli, including those in 
the RR sector and the Wilson lines. 

We derive in section |2] the free energy of a gas of weakly coupled perturbative states in 
type I superstring, in the simple case where the internal space is a factorized torus. We 
describe the induced cosmological evolution and find the radii moduli are flat directions of 
the thermal potential. In section [3| using the dual heterotic model at weak coupling, we 
correct this naive analysis by taking into account contributions of non-perturbative states 
to the free energy. In particular, D-strings modes are found to be light when the radii are 
in a neighborhood of V^") where Ai ^ 1 is the ten dimensional type I string coupling. They 
produce local minima of the thermal potential which are responsible for the stabilization of 
the radii at Iii type I? this dynamical effect occurs at strong (weak) coupling when 

D > 7 {D < Q). However, since the BPS masses of the light Dl-branes are protected by 
supersymmetry, our results are also valid at small string coupling for D > 7. In section |4| we 



reexamine the form of the corrections to the free energy along the lines of 24 , and interpret 
the non-perturbative contributions as arising from "thermal El-instantons" . What is meant 
by this is that the Euclidean worldsheets of the Dl-branes wrap the Euclidean time circle. 
In section [5} we generalize our results : The one-loop heterotic free energy is computed, 
with all of the internal moduli taken into account. We find that at certain points in moduli 
space, all scalars, except the dilaton, may be stabilized for D > 4|^ On the dual type I side, 
the non-perturbative effects induce a stabilization of the internal NS-NS and RR moduli 



^Additionally, for D > 5 the dilaton approaches a constant finite value at late times and the cosmological 
evolution is radiation dominated. For = 4, the dilaton decreases logarithmically with cosmological time 
and the coherent motion of all moduli is such that the metric evolution is that of a radiation dominated 
universe, oc However, non-perturbative effects from NS5 or D5-branes in the heterotic or type I 

theories should be taken into account in four dimensions and may play a role in stabilizing the dilaton. 
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in the closed string sector, and the Wilson lines in the open string sector. For the special 
case of D = 6, the internal volume modulus remains a flat direction, while the dilaton is 
stabilized at a small value. In section |6j we give explicit examples of loci in moduli space 
where only the fiat direction of the dilaton survives. Section [7] is devoted to our conclusions 
and perspectives. 



2 Naive perturbative type I thermal cosmology 

In this section, we derive the cosmology induced by thermal effects in the purely perturbative 
type I superstring theory toroidally compactified down to D > 3 dimensions. We shall see 
in the next section how light solitonic states correct this picture in a drastic way. In the 
following, quantities are denoted in the type I context with subscripts I and, throughout 
this paper, "hatted" ( "un-hatted" ) ones are referring to the string (Einstein) frame. Finite 
temperature Tj is implemented by considering an Euclidean time of period Pi = 2'kRiq = 
1/Ti, and coupling the associated S^{Riq) lattice of zero modes to the fermion number. We 
restrict for the moment our study to the case of a factorized internal space Y[1=d ^^i^^i) 
and analyze the dynamics of the radii Rn. 

Working in a perturbative regime, there are four contributions to the Euclidean one-loop 
partition function needed to express the free energy density, namely the torus, Klein-bottle, 
annulus and Mobius strip vacuum-to- vacuum amplitudes T, /C, A and A^. In units where 
a' = 1, a little work yields (see the appendix). 

Pi L A>0,4>0, m,n ^ i=D ^* ^ j 

A — A — rh-n 
(A,m,n)#(0,0,0) 

where Vj is the regularized volume of the {D — l)-dimensional space, cd is Stefan's constant 
for radiation in dimension D and the function G is defined in terms of a modified Bessel 



function of the second kind, Kd (x) : 

2 



ko ko 

The integer sa {sa) counts the degeneracy at oscillator level A (A) on the left (right) -moving 
side of the worldsheet, while rrii (n*) labels the momentum (winding) number along the i-th 
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cycle of the internal toru^ In (2.1), the first term in the braces is the contribution of the 



massless modes, with quantum numbers {A, m, n) = (0, 0, 0) and associated to the A/io = 1 
supergravity multiplet in ten dimensions. The Klein-bottle contribution JC vanishes. The 
annulus plus Mobius amplitude takes in a similar way the form 




Qp \ 2 — ^ 2 V 

\ A>0,m ^ 

(A,7a)5^(0,0) 

(2.3) 

where = 32 and the first term is associated to the A/io = 1 S'0(32) super- vector multiplet 
in ten dimensions. The partition function is given by the sum Zi = T + /C + ^ + A^. At high 
temperatures, it becomes ill-defined. Examining T, one finds that winding modes along the 
Euclidean time circle become tachyonic when _Rio < -Rih, where Rm = is the Hagedorn 
radius. This divergence of Zi is not a sickness of the theory, but rather the signal of a phase 
transition ||25i]. From now on, we restrict ourselves to temperatures below Tih = l/(27ri?iH). 

The free energy density is defined in terms of the partition function as = —Zij (/3iVi). 
It is expressed in terms of the G-function, whose arguments are the ratios of the spectrum 
masses to the temperature. Since 

G{x) = cd- ^ + 0{x^) when a; ~ , G{x) ~ (^] ' e"^ when a; > 1, (2.4) 
47r KzTi / 

the dominant contribution at low temperature (compared to the string scale) arises from the 



first terms of (2.1) and (2.3) and corresponds to the free energy density of thermal radiation. 



fi = - (^sl + ^ so^ CDff + ---. (2.5) 

However, if some Ru is large (small) enough, Rn > 2n Riq {Ru < l/(27ri?io)), pure Kaluza- 
Klein (winding) modes yield additional terms of the same order. The contributions associated 
to the remaining states are exponentially suppressed. 



It is straightforward to apply the techniques introduced in 15 , 16 for closed strings to 
study the backreaction of the type I free energy on the originally static background. For 
arbitrary initial conditions at the exit of the Hagedorn era, one finds that the system is 
attracted to a radiation dominated evolution, where all internal radii and the dilaton are 



*Note that the condition A — A = rn ■ n provides the level matching. 
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frozen at non-specific values depending on the initial data. Quantitatively, the final constant 
values of the R^s sit in the range 

1 



< Rii < 2'kR 



10 



D 9, 



(2.6) 



where Riq is increasing with time, corresponding to an expanding and cooling universe. 
Actually, if at some time t a radius Rij is outside this range, we find Rij{t) and Rio{t) always 



evolve so that the condition (2.6) is finally satisfied, after which the evolution of Rij comes 
to a halt. This may be seen by examining the force on the modulus fij = ln(27ri?io/-Rij) (or 



ln(27ri?ioi?y)) [15,16 



A difference compared to the type II and heterotic string cases, is that the open string 
sector is not invariant under T-duality, -Rij — )■ l/Rn (for any i), due to a lack of winding 
quantum numbers in the open sector. For instance, for arbitrary Rij (for a given j), while 



the other radii satisfy (2.6), the effective potential for Rij, which is exactly the free energy 
density, simplifies to 



= -T{ 




So CD 




+ 0{e 



2itRio < Rij 



Sq Cd + Sq 



cd + ^ G{2TxRiQ\nj\Rij 




and is shown in figure [T| in Einstein frame. When i?y < 1, the theory is actually better 
understood in the T-dual type I' picture obtained by inverting Rij. More importantly, there 
is no local minimum of the free energy density where Rij (as well as the -Rk's) can be 
attracted and stabilized. This is contrary to the heterotic case, where enhanced symmetry 
points exist and imply a local increase of the number of massless states. However, we shall 
find that the above purely perturbative analysis is missing important contributions from 
massless solitons. 
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Figure 1: Thermal effective potential (in Einstein frame) for Rij, when all other internal radii 
satisfy (2.6). The dashed curve takes only into account the perturbative type I states. The solid 



one is obtained by heterotic/type I S-duality and receives corrections from light D-string modes. 

3 Heterotic/type I cosmo logical duality 

Given that heterotic and type I theories at zero temperature are S-dual in ten dimensions, 
it is a simple but non-trivial fact that they remain S-dual at finite temperature. Technically, 
the backgrounds used to analyze the thermal ensembles are freely acting orbifolds, obtained 
by modding out with {—l)^6o, where 6o is an order-two shift along the Euclidean time circle 
and F is the fermion number. Using the "adiabatic argument" of [6] , after such a free action, 
the two theories remain dual. Since the cosmological evolutions we study are quasi-static, it 
is valid to apply at each time an S-duality transformation on the heterotic side, in order to 
derive non-perturbative contributions to the type I free energy and its resulting backreaction. 

S-dual SO {32) heterotic string 

Let us apply this point of view to the type I background considered in section [2j The 
dual theory is the SO{32) heterotic string compactified on HLd '^'^(-Rhi), where we use the 
subscript h to denote heterotic quantities. As in the type I case, the partition function 
is only well defined when the temperature Th = l//3h = l/(27ri?ho) is below the heterotic 
Hagedorn temperature, i.e. i?ho > RhH = (1 + V^)/V^- As shown in the appendix, the 
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heterotic partition function can be brought into a form divided in three parts as follows : 

9 



= I Sobo CD + y^ 2so6-i G ( 27ri?ho -Rhi I ) 

+ E SAbAG(2nR^o[^A + f2{^,-n'R^jyy)], ^^'^^ 



i=D 

9 

V L ^ ^-Tl-hi 

A>0,A>-l,m,n ^ j = D ■' 

A — A=irt-ri 
{A,m,Ti)^(0,€e'^,eei), 
Vi,Ve = -l,0,l 

where the degeneracy 6^ of the right-moving bosonic string oscillator modes is defined 
from level —1. The first contribution in Zh is associated to the massless states labeled 
by (A, m,n) = (0,0,0). They arise from the A/io = 1 supergravity and S'0(32) super-vector 
multiplets in ten dimensions. The second contribution comes from modes whose masses 
can vanish at particular values of the internal radii. For each z, these states are labeled 
as (A, m,n) = (0,eej,eej), where e = ±1 and Cj is the unit vector in the direction i. 
They are massless at the self-dual point = 1, where they enhance the gauge symmetry. 



f/(l) — 7- SU{2). The last line in (3.1) arises from the states which are never massless. It 
becomes substantial when Kaluza-Klein (winding) states become light, in the regime where 
some -Rhi's are large (small) compared to 27ri?ho (l/(27r-Rho))- AH other modes, being al- 
ways super heavy as compared to the temperature scale, yield to exponentially suppressed 



contributions (see Eq. (2.4)). 
Duality map 

In ten dimensions, the heterotic/type I S-duality identifies the Einstein frame metrics and 

,(10) n 

inverts the string couplings. Ah = e'''^ = e = 1/Ai In lower dimension these 
relations translate into a dictionary between the Einstein frame metrics, the internal radii 
and dilatons [2] : 

Rm = -7= = Rii -ju ' I = or D, ...,9, 

(n;=D2»fi«) (3.2) 

i=D 

where the Z)-dimensional dilatons are defined as (f)\^j^ = — | ^^^^ ln(27ri?h,ij)- Note 
that the Euclidean radii Rio and i?ho are included in the above relations. The inverse maps. 
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which relate the type I fields in terms of heterotic quantities, are obtained by exchanging 
the subscripts h ^ I. 

We consider non-trivial evolutions for the Einstein frame metric, dilaton and internal radii 
moduli. It is easily checked that the tree level heterotic and type I actions match, under the 
S-duality transformation (3.2) {i.e. S^^'^^ = S^^'^'^). The one-loop finite temperature effective 



potentials were computed using Euclidean backgrounds with laps functions -Rh,io in the string 
frames. For the Lorentzian Einstein frame metric, 



(3 



(27r)2 

2_ .(D) 

27ri?hioe 0-2^,1 



D-l 



y _ 2(P-1) 



(3.3) 



the corresponding first order correction to the tree level action S^^f^ is given by 



Ql-ll 



1-loop 



Jd-xV^^., Where 



(3.4) 



Note that we do not distinguish between the heterotic and type I inverse temperature /3 



and scale factor a in (3.3), as they are measured in Einstein frame and are identified under 



the duality map (3.2). To be exactly equivalent, the effective actions should be corrected 



to all orders in perturbation theory and include non-perturbative effects as well. In the 



AD) 



following, we will consider the heterotic point of view at weak coupling, e^^^ <^ 1, restrict 
our computations at the one-loop order, and deduce the type I behavior in the dual regime. 

Dual type I cosmological evolution 



To start, we apply the duality map to (3.1 ) and note that the first term exactly matches the 
sum of the first contributions in (2.1) and (2.3). This follows from the equality Bq = sq + 
{N^ — N)/2 and is due to the fact that the supergravity and 5*0(32) super-vector multiplets 
are short, with protected vanishing masses. Next, we concentrate on the interpretation and 



cosmological implications of the remaining terms in (3.1). In the weakly coupled heterotic 



AO) 



string, 6*^11 <^ 1^ the time evolution of the universe for arbitrary initial conditions at the 



exit of the Hagedorn era can be analyzed along the lines of Refs 15 , 16 . We first summarize 



the results here and then derive the dual type I cosmological behavior : 

• When all radii satisfy \Rhi — 1/-Rhi| < l/(2vri?ho)) i = -D, 9, the heterotic free energy 
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density derived from (3.1) takes the form : 



I i=D ^ 



Rhi 



-R 



hi 



C(e 




R 



h9 



1. The 



Thanks to the properties (2.4), the states with quantum numbers (A,m,n 
are responsible for the existence of a local minimum of at RhD — ■ ■ ■ — 
internal radii can be attracted and stabilized at this SU{2Y^~^ enhanced symmetry point. 
Moreover, for > 5 the string coupling e'^'t^ ' (and thus Ah) freezes to some constant value 
e'Pho determined by the initial conditions. For D = 4, the dilaton does not converge 
to a constant but instead decreases logarithmically with cosmological time. We show this 
in section |5] in a general context where we take into account all internal moduli. The rest 
of this section is valid for D > 5, while for D = 4 one has to keep in mind the late time 
evolution of 



Applying the duality map (3.2), the ratios of the masses of the above winding-momentum 
states to the temperature become : 

Mm 



27iR 



hO 







R\i 


1 


Rhi - 
j^hi 


= 27ri?io 






a7 


Rii 



(3.6) 



From the type I point of view, the corresponding BPS states have a natural interpretation as 
D (or anti-D)-strings wrapped once along the circles S^{Rii), with one unit of momentum. 
The heterotic cosmology translates into the type I context as follows. Whenever the type I 
radii start out in the dual range {Ru/Xi — 1/-Rii| < l/(27ri?io), the light D-string modes can 
stabilize them at the point 



Rj 



9 



(3.7) 



where Aio = 1/Aho ^ 1 is the late time constant value of the string coupling in ten dimen- 
sions. This implies the open string cosmology is well understood in type I, rather than in 
the T-dual picture in type I'. At each time, the width of the symmetric well of the potential 
for In Rii is ^/Xl/ {2ttRio) (see figure In total, if we denote by ^j^'' the asymptotic value of 



the type I dilaton in D dimensions and use the inverse relations (3.2), the moduli are found 
to converge as follows. 



\ (10-g)(O-2) ' 
[ZTfj 8 



Riiit) 



I AD) 10-D 



1 



e2<Pho (27r) 4 



(3.8) 
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while the temperature and scale factor asymptotic behaviors are those of a radiation dom- 
inated era, T~^{t) ~ a{t) ~ t^/^, where t is the cosmological time. Some remarks are in 
order : 



o For D > 6, (3.8) shows that the type I cosmology is at strong coupling. In this regime, 
solitons are generically light and the need to include their effects in the low energy 
effective action is not surprising. 

o For D = 6, the asymptotic values of the moduli are 6*^10' = l/(27r)^ and Rii{t) — )■ 
e~i'^M /(27r). The type I picture is perturbative. 

o For D < 6, the type I cosmological evolution is at weak coupling. However, we observe 
the necessity to take into account the contributions arising from solitons which are 
light, when we sit in the neighborhood of the enhanced symmetry points. 

In summary, for D 7^ 6 in type I, the internal radii are stabilized while the dilaton 
freezes somewhere along its flat direction. On the contrary, for D = 6, the dilaton is 
stabilized, all complex structures Rn/Rij are also stabilized at one, while the internal space 
volume n^=_D(2^-^i«) freezes along a fiat direction. This is not a surprise, since in D = 6 the 
heterotic/type I duality exchanges internal volumes and string couplings : nLz)(2^-^h,ii) ^ 



(6) 



• If at some epoch one of the heterotic internal radii satisfies i?hj > 27ri?hO; while the 
9 — D remaining ones are stabilized, i?hj = ^ ioi i ^ j, the free energy density deduced from 



(|3.1|) becomes 

CD+ V G(27rRj-^ 



J-h = -T^(so6o + (9-/^)2so6_i) 



hi 



+ 0(e-2-«'^o). (3.9) 



We see that in addition to the massless supergravity and 5*0 (32) super- vector multiplets, 
there are also contributions coming from their Kaluza-Klein descendants, which are light 
since i?hj is large. Applying the duality rules and comparing to the perturbative type I result 



in the first line of (2.7), we observe a match up to an additional contribution {9 — D) 2so&_i to 
the overall numerical coefficient. This discrepancy arises from the extra massless D (or anti- 
D)-strings responsible for the stabilization of the -Rk's at \/Xi. Therefore, the main difference 
with the pure perturbative analysis is that the plateau of the effective potential is lowered 



12 



and that the slope for Rij > 2nRio is steeper (see figure [T]). The cosmological evolution is 



however similar to the one discussed below (2.6). As their heterotic counterparts 15,16 



Rij(t) and Rio(t) evolve such that the regime where R\j{t) < 2tt Ri^it) is reached. After that, 
Rij freezes along its plateau or is stabilized at y/Xi as explained before. 

• In a similar way, if a heterotic radius satisfies i?hj < l/(27ri?ho)) while the others are 
stabilized at their self-dual points, Ru = 1 for i 7^ j, we have 



J-h = -T^(so6o + (9-/^)2so&_i) 



Cfl + ^ G (2n Ri,Q\nj\Rhj 



+ 0{e 



-27r_R,ho^ 



(3.10) 



In this case, substantial contributions arise from the winding modes along S^{Riij), which 
are light since i?hj is small enough. Their effect is to attract Rhj{t) to values larger than 



l/(27ri?ho(^)) (15,16 . Applying the S-duality rules to translate this statement in the type I 



context, we find that if Rij < Ai/(27ri?ho) at some time, the evolution of these moduli implies 
we end in a regime where Ai/(27ri?ho) < Rijy after which the internal modulus freezes or is 



stabilized at yAi. Noting that the argument of the G-function in (|3.10|) becomes 



27ri?ho l^'^l-Rhi = 2nRm \n- 



LIO 



IJ I ± 

Ai 



Ml 



(3.11) 



we conclude that the above mechanism is due to two sets of towers of D-string winding 
modes along S^{Rij). The first one contains "solitonic descendants" of the perturbative 
massless supergravity and 5*0(32) super-vector multiplets. The second set is associated to 
the descendants of the D (or anti-D)-strings responsible for the stabilization of the (9 — D) 
internal radii Ru at V^- The net result of these non-perturbative light states is to render the 
type I free energy explicitly invariant under the "non-perturbative T-duality" Rij — )■ Xi/ Rij 
(see figure 



Comments 

To conclude this section, we would like to make some remarks. We first observe that under 



the duality map (3.2), the Hagedorn radii do not match. We thus infer from the perturbative 



^ Since at late times Xi{t) — > Ajp and Rio{t) — > +00, the left-boundary of the plateau of the effective 
potential of In Rij ends by being negative. This means that Rij may freeze at some value below one. In such 
a case, a T-duality Rij — >■ ^/Rij to a type V description is more suitable. In general, a mixed type I / type 
r theory may be obtained, in order to keep all internal radii larger than one. 
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heterotic side a new value of the Hagedorn radius in type I, when Ai is large : 

f %/2 for Ai < 1 

^'""-{Vk'-^ forAi»l • ^^-^2) 

From a cosmological point of view, Rm in the regime Ai(t) ^ 1 is thus a time-dependent 
scale. Note that this non-perturbative expression for Rm obtained once D-strings are taken 

,(D) 

into account can be relevant even at weak coupling, e"**! ^1. This is for instance the case 
for D < 6, when and the -R^'s reach the asymptotic value -y/Ay^ ^ 1. 

For D > 7, the stabilization of the internal type I radii at -^/Aio ^ 1 occurs at strong 
coupling, e'^'i 3> 1. However, the D-string states responsible for this effect are BPS, so 
that their masses are protected by supersymmetry. Thus, these modes remain massless for 
arbitrary Ai, when Ru = \/Xi. It follows that the type I free energy density can easily be 



determined when Aj ^ 1 and Ru ~ -\Ai- It is actually given by (3.5), once translated in 
terms of dual type I variables. The justification of this statement is based on the following 
facts. In this regime, the string coupling is weak, e'^^ <^ 1, and the contribution of the 
perturbative part of the spectrum is that of a perfect gas. Moreover, the contribution of 
the light solitons is of identical form, since 5't/(2)'s (gauge) symmetries transform them into 
the perturbative modes in the Cartan subalgebras. We conclude that the mechanism of 
stabilization of the internal type I radii remains valid at weak coupling e"^! ^1. Since this 
yields Ru = v'Aio ^ 1, the model is better described in the T-dual type I' picture. However, 
the dynamics in the intermediate regime e"^! ' ~ 1 for D > 7 (or e"^! ' ^ 1 for D < 6) cannot 
be inferred from these arguments. 

Finally, for D < 5, additional non-perturbative states may play a role in the cosmological 
evolution. In fact, D5-branes of the type I theory (or NS5-branes in the heterotic context) can 
wrap the internal manifold in analogy with the D-strings we have considerec^ It would be 
interesting to study their effects on the dynamics, which may lead eventually to a stabilization 
of the dilaton. 



^Note that these states may contribute even for D = 5. This is to be contrasted with 5-brane instantons 
at zero temperature, which require an internal space of six dimensions. 
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4 El-instanton corrections 



We have found that non-perturbative states contribute to the type I free energy density. 
In the hterature, corrections to the low energy effective action are often considered from 
another point of view, namely instantons and their stringy generalizations. For instance, 
El contributions to holomorphic couplings have been analyzed in supersymmetric cases by 



heterotic/type I duality 24 . In the present section, our aim is to reexamine the type I free 
energy from the point of view of El-instantons and single out the configurations responsible 
for the stabilization of internal radii. In this non-supersymmetric case, we want to predict 
the El corrections in type I from dual heterotic worldsheet instantons. For simplicity, we 
restrict our analysis to the case D = 9, where instantons wrap the Euclidean time circle 
and the direction 9. This is to be contrasted with the zero temperature case where El 
corrections would only arise for D < 8. We note that by a double Wick rotation, the results 
in this section may be interpreted as the zero temperature vacuum energy contribution of El- 
instantons wrapping an internal T^, with spontaneous supersymmetry breaking boundary 
conditions along one of the toroidal directions. In this case, the temperature scale T is 
replaced with the supersymmetry breaking scale M. 

Our starting point is the heterotic model of section |3} To help exhibit the worldsheet 
instanton structure of the one-loop amplitude Zh, we work in the Lagrangian formulation 



of the zero modes lattice associated to S'^(-Rho) x S^{Rhg) (see Eqs (A. 10) and (A.ll)). We 
consider i?h9 > 1 and parameterize the zero modes by the matrix J\4. = ^ j. The case 
i?h9 < 1 may be obtained by T-duality. We may decompose the lattice sum under orbits of 
the SL{2, Z) modular group as follows. For any set of modular covariant functions /x(t, f) 
such that /x(M(r), M(f )) = /xM(r, f ), for all M G SL{2, Z), one has 

f d^-r V-'. . -X f d^'T ^ r / -X (4-1) 

n®>m9>0 



This is easily shown by applying Eq. (A. 6) twice : First to the sum over (n°,mo) and then 
to the sum over (n^, mg). The integral over the upper half plane C+ is obtained for ra^ > 
by writing mg = kn^ + / (0 < / < — 1) and changing r — )■ r + fc. The integral over T 
corresponds to the zero orbit {i.e. = 0), while the integral over 5+ corresponds to non- 
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vanishing degenerate matrices [i.e. with det = 0). The last integral over C+ is associated 
to non-degenerate matrices. 



Applying (4.1) to the heterotic partition function Zh, the contribution of the zero orbit 



vanishes due to supersymmetry, so thalj^ 



Zh- + 



nd 



zi 



7nd 



/3hVh /" d^r r(o,i6) 



(27r)9 2r| ^^fj^^ 
/3hVh /" d^r r(o,i6) 



itg > e ^2 e ^2 

mo.mg 



(27r)9 2r| r/^f 4 



g r2 g ''"2 



K,-(-l)"^°5, 



n^>m9>0 



Performing the r- integrations, the degenerate part Z^ can be brought into the form 



(4.2) 



hi ■ 



So&O Cg + ^ SAhA G i 27ri?h0 
A>0, mg ^ 



AA 



nig 



R 



h9 



while the non-degenerate contribution Z^'^ can be written as, 



7nd 



A>0, A>-1 
n^>m9>0 



9 



G'(^27ri?ho 



4A + 



h9 



h9 



(4.3) 



(4.4) 



Summing over mg in (4.4) enforces the level matching condition A — A = n^mg for some 
integer mg, whenever ^ 0. The "missing term" for = is actually the contribution 



of the degenerate orbits Z^. In total, Z^ + Z^'^ yields with no surprise the expression (3.1), 
which can be analytically continued in the range 1 < R^g < V^- However, to exhibit the 
instantonic structure, it is better to leave the sum over rhg. 

In Z^, only pure Kaluza-Klein modes along the directions 9 and contribute and the 
worldsheet embedding in the target torus is trivial (no instanton number). Therefore, these 
states do not play a role in stabilizing the internal circle. In order to extract the configura- 
tions in Z^"^ responsible for fixing R\^g at the self-dual point, we know it is enough to focus 



^ The use of Eq. ( |4.l[ ) is valid if the argument of the discrete sum to integrate is absolutely convergent. 
In the present case, since the right-moving block r^o ig)/^^'* and the left-moving Og,/rf character involve 



diverging powers of e^'^'^'^ in the hmit T2 — >■ -l-oo, Eq. (4.1) can be trusted if i?ho > a/S and i?h9 > V2- The 



first condition is not problematic as we are focussing on the dynamics at low temperature. Since we are 
interested in the stabilization of i?h9 around 1, the second condition could be a problem. However, we see 
shortly that the final expression (4.4 1 can be analytically continued all the way to i?h9 = 1- 
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on the dominant contributions in the low temperature expansion. The terms with A > 1 
are exponentially suppressed, 0{e~^'^^'^'^), compared to the contribution with A = 0. The 



latter arises from BPS configurations and, at this level of approximation, Z^'^ in Eq. (4.2) 
involves a purely antiholomorphic function, B{f) = r(o,i6)/?/^^, dressed by an inverse power 
of T2 and the lattice of zero modes associated to the directions and 9. This form is similar 
to the one encountered in the evaluation of holomorphic couplings, when supersymmetry is 



unbroken 24 



We can now define instanton configurations, with associated Kahler and complex struc- 
ture moduli T and 3^ as, 

T = tT2= i{2ko + l)Rm ■ n'Rh 

y = yi + 



Instanton with > rrig > , A;o > : 



mg , {2ko + l)Rho 
+ 1 



(4.5) 



n9 n^Rh9 

where (2fco + 1) is the instanton number, which counts the number of times the worldsheet 
wraps around the target torus. Using these notations and introducing coefficients q;„ G N in 



the expansion of the Bessel functioE^ in (2.2), K9{x) = A/7r/(2a;)e ^^n=o^n/x"', we may 
write (4.4) as 



7nd 



(10) 



{2tt 



So 

instantons 



,10 r^yi 



n=0 



(27rT 



A>-1 



4— n 



+c.c.+0{e' 



(4.6) 

where V^^^^^ is the ten-dimensional Euclidean volume. This result can be given a more elegant 
appearance by noting that B{y) is a modular form of weight 4. Introducing the modular 
covariant derivative DX = {dy + where X{y) is any modular form of weig ht 

the brackets in (4.6) become l/(7rT2)" X]m=o 7""i(^^2)'^-D™'S(3^), where 7„m are rational 
numbers. 

The above expression of contains far too many explicit terms needed to study the 
stabilization of i?h9- In (4.4), the dominant contribution for A = arises when A = —1 and 

= 1, while the remaining terms are exponentially suppressed, 0{e~'^^^^°). Restricting to 
A = —1 and the instanton configurations = 1, rhg = 0, fco > in Z^'^, we can add the 
degenerate contribution = (/3hVh//3^) so&oCg + 0{e~'^'"^^°) to recover the first line of Eq. 



(3.1) required for the derivation of the stabilization of R^g. 



^In any odd dimension, the Bessel function admits a power series with a finite number of terms. 
^This means that X{y + 1) = X{y) and X{-l/y) = Xiy)/^^. Moreover, DX is a modular form of 
weight r — 2. 
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We now wish to interpret Eq. (4.6) from the perspective of the type I superstring. Under 



the heterotic/type I dictionary (3.2), the complex and Kahler structures 3^ and T are mapped 



into 3^1 and Tj/Aj. Consequently, the exponential factor of T in (4.6) yields the exponential 
of the Nambu-Goto action for a D-string, and Z^'^ translates into a sum of El instantons as 



m 



24 



El 



(10) 



e 



So 



El instantons 



12 n=0 



(27rT 



12 



A>-1 



A—n 



+ c.c. + 0{e 



-An 



Rio 



(4.7) 

Actually, the configurations of the D-string worldsheets wrapped on S'^(-Rio) x S'^(-Ri9) are 
highly dissymmetric at late times in the sense that Rio{t) — )■ +oo and Rigit) ~ a/ Xi{t) — t- 
a/Aio- However, this does not mean it is artificial to consider such El-instantons. Instead, 
they open the possibility to derive from a pure type I point of view the free energy responsible 
for the stabilization of the internal moduli (or the effective potential at zero temperature 
when at least two internal directions are compactified and supersymmetry is spontaneously 
broken). Thus, it would be interesting to derive D-brane instanton corrections from first 
principles, in the case where supersymmetry is spontaneously broken. The full instantonic 



structure of (4.4) should also be interpreted from a type I point of view, even when all 



contributions with A > and A> —1 are kept explicitly. 



5 Heterotic and dual type I moduli stabilization 

We would like to extend the analysis used in section |3] to include the remaining moduli in 
addition to the internal radii. We consider the heterotic string compactified on T^^~^ at a 
generic point in moduli space and show that when finite temperature is switched on, the 
free energy density can stabilize all internal moduli. Our study is based on the effective 
action at finite temperature and weak coupling for the massless degrees of freedom, while 
all massive states are integrated out. Introducing simplified notations, we are interested 
in non-trivial backgrounds for the Einstein frame metric g, the dilaton (p in D dimensions 
and all real- valued internal moduli which we denote collectively as $. Concretely, $ 
contains the components of the metric Qij and antisymmetric tensor Bij, together with the 
Wilson lines Y/ {i, j = D, . . . ,9; I = 10, 11, ... , 25). It is then straightforward to deduce the 
dynamics and final expectation values of the type I counterparts of these scalars by using 
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the duality map 

9ij = ^ J = Cij , F/ = Y[{ , (5.1) 

M 

where Cij is the RR 2-form. Detailed examples of this analysis will be given in section [6] for 
D = 8. 

The heterotic low energy effective action 



S = I d'^x^ 



(5.2) 



involves the tree level moduli space metric Fmn = Fmn{^) and the one-loop free energy 
density J-" = J^{T, 0, $). Since the backreaction of J-" on the classical background is already 
a one-loop effect, there is no need to take into account the quantum corrections to the 
kinetic terms. For homogeneous and isotropic evolutions, variation of S with respect to the 



time-dependent metric (3.3), dilaton and moduli yields, in cosmological time defined by 

dt = /3{x^)dx°, 

~ ^^^^ ~ H' = ^^^' + \ Fmn^''^'' + P, (5.3) 
~ ^^^^ ~ H' + {D- 2)H + 02 + \ Fmn^''^'' + P = 0, (5.4) 

+ (D - l)H^ + 7^ = 0, (5.5) 

<|.^ + (D - 1)//$^ + F^-^^ (f^pq - ^FpQ^) <i.^<i.« + F'^'^Tj, = 0. (5.6) 

In these equations, H = d/a and the thermal pressure and energy density are found to be 

dP 

P = -J', p = T—-P. (5.7) 

Additional indices and denote partial derivatives with respect to and while 
jpMN ^ ^jp-i^MN ^ convenient to replace Eq. (5.4) by the constant entropy constraint. 



The latter is found by integrating the energy-momentum tensor conservation law derived 
from the above differential system (see [15]), 



p + (L>-l)iJ(p + P) = 0J'0 + $^J'm =^ a^"^ = constant entropy. (5.8) 

In order to find particular evolutions characterized by static moduli, (0(t),$(t)) = 
bQ,^o), we need to specify J-". For any supersymmetric spectrum, the one-loop free en- 
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ergy density is 

^ - -e-^ n ^ E E (1 - i-r^) --t^Hg f ^ 

Jo 2/ (27r/)f ^ K ^ ) ) ^ \ T 

(5.9) 

where is the mass of each boson/fermion pair s, and the dilaton dressing in front of 
the integral is introduced to switch from string to Einstein frame. This general expression 
applied to our case of interest, namely the heterotic string on T^^^^ , is explicitly derived 
from a one-loop vacuum-to-vacuum amplitude in the appendix. In the notations introduced 
there, soTq = 2^ x 24 boson/fermion pairs of states are massless everywhere in moduli space 
while the other modes have moduli-dependent masses, Ms{^). Since light states have the 
tendency to lower J-", effective potential wells can be found at any point $o where uq > 
pairs of modes generically massive are accidentally massless, M„($o) = 0, u = l,...,nQ. 
The fact that we have at zero temperature 16 real conserved supercharges implies that such 
points are associated to enhancements of the gauge symmetry. Defining Mmm to be the 
lightest non- vanishing mass at $0) the free energy density can be written in a neighborhood 
of $0 as, 

^^-r" (5.10, 

At low enough temperature, the exponentially suppressed terms can be neglected and we 
may derive identities for the thermal source terms at $05 including the equation of state, 

^1^^ = -T^(soro+no) CD , J-^l^^ = , TmU^ = , p|^„ = {D-1)P\^^ oc T^. (5.11) 

It is then straightforward to check that the evolutions 

ao(t) (X (X t^/^ , <Pit) = <Po, m = %, (5.12) 

corresponding to radiation eras with static moduli are particular solutions of the equations 
of motion. 

The above trajectories are actually attractors of the dynamics in some circumstances. 
To study this, we analyze their stability under small time-dependent deviations, 

a = ao(l + e,), T = To(l + er) , = </)o + , $^ = $f + e^. (5.13) 
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They are associated to the supergravity and super-vector multiplets of the 5*0(32) Cartan generators. 
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We first perturb the internal moduli equation (5.6). Denoting Hq = ao/ao, one obtains at 
lowest order, 



■■M 



{D-l)Hoi 



M 



iV AT e 



where 



N 



:^ML\ 



LN 



{To,0o,*o)- 



(5.14) 



A V is an effective "time-dependant squared mass matrix" evaluated for the background 



(5.12). Since 



T 



D-2 



A4l_ Cd-2 

gD-2 

4:71 



no 

E 

u=l 



N 



(5.15) 



$0 



is semi-definite positive, A V is diagonalizable with non-negative eigenvalues , which we 



define as 



£)2f2(D-2)/D ■ 



In the case when some Am's vanish, one needs to take into account 



quadratic terms in Eq. (5.14) (see the discussion of the dilaton equation below). In par- 
ticular, this is required when moduli sit on the plateau of their thermal effective potential 
(see figure [T]). For simplicity, we proceed by analyzing the most interesting case, where 
all internal moduli are "massive", which means A^^ > 0. Switching to a diagonal basis of 



perturbations e , one obtains from (5.14) 



(5.16) 



where are integration constants and J:d^ are Bessel functions of the first kinc 



This 



describes damped oscillations with amplitude of order l/yt, where t is supposed to be large 
enough so that |e^^| ^ 1 is satisfied. 



Next, we derive from (5.5) the equation for the dilaton perturbation at leading order, 

4 



where J^, 



<f>MN I $Q 



D~2 



J^MN\$g- (5.17) 



Since the constants C^^ are a priori of order one, we take into account the quadratic source 
in "massive" epsilons. Thus, can be written as the sum of the general solution to its 



homogeneous equation, plus a particular solution to Eq. (5.17). The former is of order 
l/aQ~^ and turns out to be dominated at late times by the latter. Actually, using (5.16), 



^^This follows from the fact that the matrices F~^^^ and T are (semi-)definite positive, so that 
F^^/^^FF^'^/^ ~ ^1/2^^-1/2 jg semi-definite positive. Note that in models where the spontaneous breaking 
of supersymmetry is generic i.e. not only due to thermal effects, each term in the sum over the boson- fermion 
pair u in Eq. (5.151 is dressed with a -|- (or — ) sign when the boson (fermion) is lighter than the fermion 
(boson). In such cases, F is not semi-definite positive and the extrema of F can be minima, maxima or 
saddle points. 

-'^^For D — J_i should be replaced by the Bessel function of the second kind, y_i. 



21 



the quadratic source term involves products of Bessel functions with arguments Xpt^^^ 
and Agt^/^. Integrating it once, the dominant contribution to a^'^e^ is found to arise for 
"constructive interferences", i.e. when Ap = \q. This yields the asymptotic behavior, 

Cs 1 



,D-2 



e<t> oc 



^l-4/D 



for D > 5 and oc Int for D = A, 



(5.18) 



where is a fully determined coefficient quadratic in C^^'s and positive. For D > 5, the 
consistency condition \e^\ <^ 1 is automatically fulfilled at late times. On the contrary, 
the case D = A yields formally to a logarithmically decreasing and one may worry that 
the our expansions breaks down. Therefore, we have solved numerically the full non-linear 



differential system (5.3)-(5.6) in this case and found that the perturbative analysis gives the 
correct late time behavior, which we summarize at the end of this section. 

To analyze the evolution of the scale factor and temperature fiuctuations, we expand 



the energy density and pressure around the background (5.12) and find from Friedmann's 



equation (5.3) and (5.8), 



1 



D-2 



{D -1){D- 2) Hoea=- FMNk/''e'' + Dp\ 



D-3 



(To,</.o,*o) 



MN 



(To,</.o,*o)^ ^ ' 



(5.19) 
(5.20) 



It is then straightforward to solve for e^, whose asymptotic behavior is again dictated by 
the source terms in "constructive interferences" arising from the products e^^e^ and e^^'^e^ 



in (5.19) and (5.20). The late time scaling property of is found to be 

a| 1 



(5.21) 



which can be used in Eq. (5.20) to find 



er oc (1 + oscillations with constant amplitude). 



(5.22) 



In (5.21) and (5.22), the coefficients of proportionality are again fully expressed in terms of 
the s. 

We signal that for D > 5, all terms we have neglected in the perturbed equations of 
motion are a posteriori found to be dominated by the sources we took into account. This 
guarantees the validity of the asymptotic behaviors we have found for the deviations defined 
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in (5.13). These results have been confirmed by direct numerical analysis of the unperturbed 
system of differential equations in some examples. Since all fluctuations converge to zero, 
the late time cosmology is radiation dominated. In particular, the dilaton motion and the 
damped oscillations of e^^ store a negligible amount of energy as compared to the thermal 
radiation energy. The internal moduli are dynamically stabilized and their effective time- 
dependent masses (measured in Einstein frame) are M^m oc Tq ^ . 

For D = 4, the numerical simulations show that the internal moduli converge to $05 
while the dilaton decreases logarithmically with time. Individually, the energy stored in the 
dilaton motion, the total energy (kinetic plus potential) of the damped oscillations of $, and 
the thermal radiation energy decay at the same rate. Their late time behavior satisfies 

so that the metric evolution is identical to that of a radiation dominated universe, a oc 

The above logarithmic behavior of the heterotic dilaton is transferred by heterotic/type 
I duality to the type I dilaton for D = 4. Moreover, in any dimension, stabilization of the 
internal moduli on the heterotic side implies stabilization of internal moduli on the type I 
side, except for the special case of D = Q, where S-duality exchanges the six- dimensional 
heterotic coupling with the type I internal volume modulus. 

6 Example : Dual heterotic/type I strings on T"^ 

Our aim is to illustrate the analysis of the previous section with examples for D = 8. We 



want to find local attractor solutions of the form (5.12) associated to enhanced symmetry 
points $0 of the internal moduli space of the heterotic string on T^. We shall see that the 



one-loop free energy density has enough structure to stabilize T = Bgg + i\/ g&sQm — S'fg, 
U = {g89 + i V 988999 - gig) / §88 and the Wilson lines Y/ = 8, 9; / = 10, 11, ... ,25). This 
translates in the type I side into expectation values of the closed and open string internal 
moduli via the duality map T = Ti, U = Ui, Y/ = where 

^ , . V 91889199 - 9189 , ■ - 6j {91889199 - gimf^^ -.. 9\89 + «a/ 9\889\99 " .^189 
/I = (--89 + * ^ = <-'89+^e , U\ = 2 • 

Ai Zn gis8 

(6.1) 
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The only remaining flat direction of the thermal effective potential corresponds to the 
heterotic and type I dilatons in eight dimensions, which are related as : (pi = —^(p — 



fin ((27r)V^88^99 -^Ig)- 

The heterotic effective action in the Einstein frame is (see for instance appendices D and 
E in (26)) 



S 



R [dcPf I (\dUW\dT + Y{^dY^^^^\udYi -dY^\''^ 



4 \ W 



V 



T2U2 



7\. 



(6.2) 



Indeed, if we arrange the thirty-four entries of the moduli vector as $ = {Ti,T2,UiM2,Y^ , V/'), 
where indices 1 and 2 refer to real and imaginary parts, the metric components of the general 



expression (5.2) are 



2Tf 





1 

2Ti 



sym. 









yj 



2Ui 



2W 









2T2U2 



\ 



\ 





2T2U2 " sr? 

1 Ki'j' ^ yi'yi' , 



2T2U2 



+ 



(6.3) 



The free energy density J-" is determined by the mass spectrum (see Eq. (5.9)), which 
is specified by the left (right) -moving oscillator number A {A), the internal momenta and 
winding numbers mj, ra* (z = 8, 9), and the root vector of the right-moving internal lattice 
rspm(32)/Z2- As reviewed in the appendix, the mass formula = 2{A + A) + ^ (p| + T^r) 
involves the left and right-moving momenta along the compact directions. 



vIr = [rn, - Q-'Y/ - B,, - ^n^Y/Y/y*'' T n'e\ for ^, j, J = 8, . . . , 9; J 



10, 



,25, 



for / = 10, . . . , 25; Q G V syin{?,2) jiL, 

(6.4) 



where = efej and e**^ej = 5*. More explicitly, one obtains 



"J 

Ml . ^JTM.Y) = 



where we have defined 



-nisU + nig + rn"" + [TU - -W'W^) + W^Q^ 



+ 4v4, 
(6.5) 



:= UYi - F/ , f:=r+ -Y^'W 



(6.6) 
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and used the level matching condition, A — A = min^ + ^Q^Q^ . At generic points in moduli 
space, the gauge symmetry is U{1)1 x f/(l)|. x U{lY^, where U{1)1 x f/(l)| arises from 
compactification, and U{iy^ is the Cartan subgroup of SO{32)r. We now examine special 
points in moduli space where no pairs of bosonic and fermionic superpartners generically 
massive are accidentally massless. Since at zero temperature the model is maximally super- 
symmetric, such points are associated to enhanced gauge symmetries. In fact, the additional 
massless modes arise at oscillator levels A = 0, A = —1, so that no is proportional to 
Sor_i = 2^ (see the appendix) and the enhancements of the gauge theory arise from the 
right-moving sector only. In the following two examples, we will simplify the notations by 
omitting the subscript "i?" in the right-moving gauge group factors. 

Local attmctor 1 : U{l)l x SU{?>) x 50(32) 

We start with the most obvious attractor where all Wilson lines vanish, = 0, leaving the 
50(32) group unbroken. The torus moduli take the values T = U = ^ + i^, implying an 
additional SU{3) gauge factor. The hq states responsible for the enhancement of ^(1)^ x 
U{lY^ SU{3) X 50(32) are divided into two groups : 

• 6 X 2^ boson/fermion pairs imply f/(l)^ — )■ SU{3). Their quantum numbers are (m, n) = 
±(1,1; 0,1), ±(0,1; -1,1) or ±(1,0; 1,0), and Q = 0. In this case, p^^ = = p^j^^^ and 
{p%,p%) realize the root vectors of SU{3), which represent a hexagon. The corresponding 6 
mass formulas are, 

' :^|l-W±rW- ^WWp for (m,n) = ±(1,1,0,1), 
Km,n,d={ i^|l-r±rW-iWW|2 for (m,rl) = ±(0,1,-1,1), (6.7) 



r2U2 



\r-U\'^ for (m,n) = ±(1,0,1,0). 



• 480 X 2^ boson/fermion pairs to recover U{lY^ — )■ 5*0(32). They have {rh,n) = (0,0), 
Q = ±(1, ±1, 0, 0), ±(1, 0, ±1, 0) or any other permutation. In this case, p^'^ = = 
p^jf , while (p/f~^°) realize the root vectors of 50(32). The corresponding 480 mass formulas 
are 

MoMQ = 7^l±("^'±"^"')r ' /,J=10,...,25, J. (6.8) 



To compute the squared mass matrix defined in Eq. (5.14), we first evaluate the second 



derivatives (5.15) of the free energy at $o- The non vanishing components are proportional 
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1, 2 (no sum over a) 



160 xr , i 



8,9; / 



10, . . . , 25 (no sum over i, I). 

(6.9) 



The nonzero entries of the metric (6.3) at $o are also found to be 

2 

3 ^ 

Fvivi = —- , z = 8, 9; / = 10, . . . , 25 (no sum over z, /). 



T T 



a 



1, 2 (no sum over a) 



(6.10) 
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The resulting matrix of squared masses is diagonal, with strictly positive eigenvalues. There- 
fore, all fiat directions of the internal moduli space are lifted. Once the dynamics is attracted 



to the trajectory (5.12), the "time-dependent moduh squared masses" are 

C6 



^23x24e^To^ 
47r ^ 



or 



M, 



$2 



2^ X 240 e^Tn^ 



(6.11) 



The first one corresponds to Ti, T2, Ui, U2, while the second is associated to the Wilson lines 
y/ and Yg . The additional factor of ten for the latter can be understood from the fact that 
they are coupled to ten times as many additional states as compared to the torus moduli. 



Local attractor 2 : U{1)1 x SU{2) x 50(34) 

The point $0 "we now consider corresponds to the values T = lA = i/v^, Y^. 



/>io 



and 



^9 



Vl2 



V25 



-1/2. This moduli configuration is much less trivial 
than the previous one, since it is going to give rise to the gauge group SU{2)^ x 5'0(34)9 25, 
where the subscripts denote which directions z = 8, 9 and / = 10, ...,25 are associated with 
the gauge factors. There are no = 546 x 2^ extra massless boson/fermion pairs of states, which 
can be divided into 2x2^ for the SU{2)^ and 544 x 2^ for the 5'0(34)9 . 25 enhancements. 
Note that the 5*0(34)9 . 25 factor arises from an enhancement of the f/(l)9 symmetry of the 
torus, with the S'0(32) symmetry of the internal lattice. The detailed quantum numbers 
of the extra states follows : 

• 2 X 2^ boson/fermion pairs give f/(l)8 — ^ 5'f/(2)8. They have (m, n) = ±(1, 0; 1, 0) and 
Q = 0. In this case, p^-® = = while = ±y/2 realize the root vectors of SU{2)s. 

For S'0(34)9_...^25, the 544x2^ pairs of bosons and fermions giving U (1)9'^ 25 ~^ 5'0(34)9^,..^25 
are subdivided into : 
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• 420 X 2^ pairs transform in the adjoint representation of 5*0(30) and are giving rise to 
f^(l)n....,25 '5O(30)ii,...,25. 210 x23 have {m,n,Q) = ±(0, 1; 0, 0; 0, 1, 1, 0, 0) or any per- 
mutation of the last 15 entries. The other 210 x 2^ have (m, n, Q) = (0, 0; 0, 0; 0, 1, —1, 0, 0) 
or any permutation of the last 15 entries. 

• 60 X 2'^ pairs transform as (2,30) under 50(2)10 x 5O(30)ii^...^255 giving the enhanced 
group SO(32)io,...,25- 30 x 2^ of them have (m,n,Q) = ±(0, 1; 0, 0; -1, 1, 0, 0) or any 
permutation of the last 15 entries. The other 30 x 2^ have (m, n, Q) = ±(0, 0; 0, 0; 1, 1, 0, 0) 
or any permutation of the last 15 entries. 

• 64 X 2^ pairs transform as (2,32) under 50(2)9 x 5*0(32)10,. ..25, giving the enhanced 
gauge group 5*0(34)9 ,,. 25- 32 x 2^ of them have {rh,n,Q) = ±(0, 1; 0, — 1; |, |) and 
±(0, 1; 0, —1; — |, — ^, ^, ^) or any permutation of the last 15 entries. The other 32x2^ have 
(m, n, Q) = ±(0, 2; 0, —1; — |, |, |) and ±(0, 2; 0, —1; — |, |, ^, ^) or any permutation of 
the last 15 entries. 



Proceeding as before, the squared mass matrix in (5.14) can be evaluated. Its diagonal- 
ization reveals two groups of eigenvalues. 

Ml = ^2^x166^ To^ Ml = ^2^x 2566^ (6.12) 

*i 47r 2 47r 

The first one is associated to Ti — Ui — \{Y^^ — Y^^ _ . . . _ Y^^) and T2 — W2, while the 
second corresponds to Ti+Ui, T2+U2 and all 32 Wilson lines. Thus, we find a second point 
in moduli space where all internal moduli are stabilized by the thermal effective potentia 
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7 Conclusions and perspectives 

In this paper, we considered toroidally compactified heterotic and type I superstrings at finite 
temperature. Applying the rules of heterotic/type I duality, we inferred novel contributions 
to the free energy of a gas of type I superstrings. These contributions are due to BPS D- 
strings wrapped on internal circles which become massless at special points in moduli space, 
enhance the gauge group, and lift flat directions. These conclusions are based on the S-dual 
heterotic picture at weak coupling. At finite temperature, the latter is a no-scale model i. e. 

^•^We have also investigated a third local attractor at the point T = U = i/2, Y^-^^ = = Fg^"'^"'^, 
y^i2,...,25 _ \i2^ which corresponds to the gauge enhancement SU{2) x SU{2) x 5*0(32) and a stabilization 
of all internal moduli. Due to its similarity, we do not present its details here. 
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a flat background where all supersymmetries are spontaneously broken at tree level. 

We computed the one-loop free energy density on the heterotic side for D > 4 and 
found points in moduli space where all internal moduli are dynamically stabilized due to the 
cosmological evolution. Additionally, in D > 5, the evolution of the dilaton asymptotes to 
a constant value, while in D = 4, the dilaton turns out to have a logarithmically decreasing 
behavior. 

Using the S-duality, this implies that for D > 7, all type I internal moduli can be 
stabilized at strong coupling. In D = 6, the S-duality maps the heterotic coupling into the 
type I volume modulus. As a result the only remaining flat direction in type I is the internal 
volume modulus, which asymptotes to a constant finite value, while the type I dilaton is 
stabilized at weak coupling. For the cases D < 5, all type I internal moduli can be stabilized 
at weak coupling. Furthermore in D = 4, the type I dilaton inherits the logarithmic behavior 
from the heterotic dilaton, while it asymptotes to a constant in higher dimensions. In all 
cases, the late time geometric evolution is identical to a radiation dominated evolution. 
Furthermore, all solutions are stable under small perturbations and are thus local attractors 
of the dynamics. 

It is worth stressing that the effects of the massless BPS non-perturbative D-strings 
persist at weak coupling, as their masses are protected by supersymmetry. As a result, the 
stabilization in type I for D > 7 persists at weak coupling. Furthermore, taking these modes 
into account is not optional in phenomenologically motivated uses of the type I superstring. 
Actually, this is not the first time massless solitons play an essential role in weakly coupled 
theories. For instance, in type IIB compactifications on Calabi-Yau threefolds, the conifold 
singularities in the vector multiplets moduli spaces are explained by massless hypermultiplets 



realized by D3-branes wrapped on vanishing 3-cycles 27 



Realistic models should include also a spontaneous breaking of A/4 = 1 supersymmetry 
at a scale M, before finite temperature T is switched on. In this case, the universe is 



attracted to a "radiation-like dominated era" 11-16]. This evolution is characterized by 
coherent motions of e^'^^^^ (where (p is the dilaton in four dimensions) and the modulus M{t), 
both proportional to T{t) such that Friedmann's equation is effectively that of a radiation 
dominated era, oc T^. The energy stored in the oscillations of the moduli around their 
minima is found to be dominated by the thermal energy and so the stabilization of the 
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scalars is guaranteed. Moreover, infrared effects are expected to put a halt to the run away 
behavior of the string couphng and supersymmetry breaking scale. In particular, when 
T{t) reaches the electroweak scale M-^w, radiative corrections are not screened anymore 
by temperature effects and the electroweak breaking is expected to take place [28] . This 
should be accompanied by the stabilization of M{t) around Mew 29 . Clearly, it is of 
utmost importance to implement these effects in our cosmological set up since this would 
provide a precise context for addressing questions of dark matter, astroparticle physics and 
phenomenology. Additionally for D = 4, as well as D = 5, there is the possibility of large 
contributions coming from light NS5-brane states in the heterotic theory or D5-brane states 
in the type I theory which have not been taken into account yet. It is possible that these 
states can play a role in stabilizing the dilaton. To make progress in this direction, one may 
try to exploit heterotic/type II duality in D = 4 which is a strong- weak duality. 
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Appendix : Thermal partition functions 

Type I superstring 

To study the canonical ensemble of a perfect gas of maximally supersymmetric open and 
closed superstrings, we compactify the type I theory on the Euclidean background ^^(i^io) x 
T^"^ X HLd S^i^ii)- Bosons (fermions) are imposed periodic (antiperiodic) boundary con- 
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ditions along S'^(-Rio), where (3i = 2tcRio is the inverse temperature. The spatial torus T^~^ 
is considered in the large volume Vi limit. Our aim is to compute the one- loop thermal par- 
tition function. The treatment of a generic Scherk-Schwarz compactification can be found 
in 



30 and the case of present interest is reviewed in 31 



In the closed string sector, the torus contribution is half that of type JIB, 
AVi I r (fr 1 



T 



^ 2t. 



§+1 r/8r/8 



(27r)^ 2 
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AVi 1 f d'^T 1 
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'2^ 



\a+b+ab ^[b]^ 



qA^LqiPR 



\fho(a+a)+rfi{b+b) 



(A.l) 



n" even, mg 



+ E 



e ^2 



■1)-0(V8S8 + S8V8: 



rfi odd, mo 

where g = e^*'^'^ and pL,ffi = rrii/ Ru =F n^Ru- The above second expression involves 5*0(8) 
affine characters, where those associated to the vectorial and spinorial representations satisfy 

V^^ ^8 
jj8 ^8 



E 

A>0 



(A.2) 



The Klein bottle amplitude /C is obtained by keeping all characters of T which are invariant 
under left -H- right symmetry. Symmetrizing and antisymmetrizing the NS-NS and RR 
sectors respectively, JC involves the combination Vs — 5*8 and is thus vanishing. In the open 
string sector, the thermal annulus and Mobius strip amplitudes are 
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where = 32, q = e 
power expansion 



Pi 



rrii/Rii and the "hatted" characters in Eq. (A.4) have the 



A>0 



(A.5) 
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We proceed by evaluating more explicitly the amplitude T by "unfolding" the fundamen- 
tal domain of integration [32]. In fact, for any set of modular covariant functions f(n,r~n){Tj ^) 



such that f{n,m){M{T), M{t)) = f(n,m)M{'T-, t) for all M G SL{2, Z), one has 
r (p^Y r (Pi r 

/ —y^f{n,rh){T,f)= —fiOfi){r,f)+ — y'/(0,m)(r,^), (A 

^2 „^ ^2 JS+ ^2 ^_^o 



6) 



where 5+ is the upper half strip : —1/2 < ri < 1/2, T2 > 0. Applied to Eq. (A.l), 



supersymmetry implies that the contribution for uq = itlq = vanishes and we are left with 



integrals over 5+ for hq = 0, rfiQ 7^ 0. Defining rfiQ = 2kQ + 1 and using (A. 2), one obtains 
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where level matching is implemented by integrating over ti. Using the formula dx- 



-a/ x — bx 



2a 2 5 2 Ky^i{2y ab)^ where K^{x) is the Bessel function of second kind, the integral over 



T2 yields Eqs. (2.1) and (2.2). Similarly, applying the expansions (A. 2) and (A. 5) in Eqs 



(A. 3) and (A.4), we have 
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which gives Eq. (2.3) after integration over T2. 



Dual heterotic string 

We proceed by deriving the partition function of the dual heterotic theory, which is com- 
pactified on S'^(-Rho) x T^^^ x HLd 'S'^(-Rhi)- Bosons and fermions are again given periodic 
and antiperiodic boundary conditions along the Euclidean time circle, whose circumference 



Eq. (A.6) is true as long as it is allowed to exchange discrete sum and integration, a fact which is 



guaranteed if the integrand is absolutely convergent. This condition is satisfied for T when i?io > i?iH- 
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defines the inverse temperature /3h = 27ri?ho- This yields 
/3hVh f d?r r(o,i6) 



ri24 



Lqi^R 



a,b 



Ml 

^4 



\^a+b+ab ^lb\ ^ •jmoa+n^b+mon'^ 



(0,16) 



qiPLq^PR 



e ^2 



lia|nOT+mo|2 



n" even, mo 



+ ^ e ^2 

odd, mo 



^^|nOr+mo|2 



;-i)'^°08 - Cg 



(A.IO) 



where g = e^*'^'^, while PL,Ri = ^i/Rhi =F 'T-'-Rhi and the volume Vh are now measured in the 
heterotic theory. Alternatively, the lattice of internal zero modes can be considered in its 
Lagrangian formulation, as needed in section |4] for the direction 9, 
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To unfold the fundamental domain of integration in (A.IO), one can use the identity (A. 6) 



as in the torus amplitude in type I. Expanding the 50(32) right-moving lattice as 
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and using Eq. (|A.2|), one obtains 
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fco , m, n 
A>a^A>-l 
A — A — m^n 



which can be integrated to give Eqs. (3.1) and (2.2) 



Heterotic string at generic point in moduli space 

In sections |5] and |6] for D = 8, we study in the context of the maximally supersymmetric 
heterotic string the stabilization of all internal moduli by the free energy density at weak 
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coupling. In Einstein frame, the latter is J-" = — e-D-2'^^j^/(/3jj\/j^)^ where is the dilaton in 
dimension D and is the vacuum energy in the Euclidean background S'^(-Rho) ^ T 



D-l 



rpio-D^ The internal moduli are the metric Qij, the antisymmetric tensor Bij and the Wilson 
lines Y/ {i, j = D, . . . ,9; I = 10, 11, ... , 25). Proceeding as before, the partition function 
takes the following forms, 
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where we introduce the coefficients of the expansion 77 = Y1a>-i^a1^- The moduli- 
dependent internal momenta are specified by m, n and the root vector of the right-moving 
lattice Tspin(32)/Z2 [33j , 

= (^i - - n^Bij - \n'Y,'Y^') e"' T n'e{ for 7, j, / = . . . , 9; J = 10, . . . , 25, 



= 72 (g^ + n^F/) for / = 10, . . . , 25; Q e Tsp^nm/z. 

(A.15) 

where {cj} is a vector basis of T^'^'^ i.e. gij = e[ej and = Sj. Since these momenta 

satisfy ^{Pl—Pr) = —2rh-n — Q-Q, the level matching condition implemented by integrating 



over Ti in Eq. (A. 14) is A — A = m- n+ \Q ■ Q, which yields 
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A>0, A>-1 
A — A—jn-n+\Q 



where ^ ^(^, 5, F) = 2(A + A) + ^ + j}fj) are the masses of the boson/fermion pairs 
of superpartners. Integrating over r2, the above expression for Zh leads to the free energy 



density (5.9), while for = 8 the mass spectrum takes the more explicit form (6.5). 
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